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On the charged meson pair theory of nuclear forces 
Part I 


By O tor Bru Lin 


_ Introduction 


The problem of giving a consistent theory for the nuclear forces has not yet 
— even in the simplest case — been finally solved. The main experimental 
facts to be accounted for are — in the simplest case of two nucleons — the 
binding energy and electric quadrupolemoment of the deuteron, the proton- 
neutron and the proton-proton scattering. The magnetic moments of the nucleons 
and the light nuclei are moreover very accurately determined and the problems 
of the saturation of the nuclear forces and f-radioactivity have to be kept 
in mind. Many of these data can phenomenologically be described by several 
different potentials; potential well, Gauss potential, meson potential, etc. 
These potentials must be spin-dependent and, at least to some extent, non- 
central. It may, however, be necessary to assume different ranges of the nuclear 
forces in the singlet and triplet state of the deuteron. To get a survey of 
the experimental facts, discussed from a theoretical point of view, we would 
refer the reader to Brtrur, Elementary Nuclear Theory [47], and ROsENFELD, 
Nuclear Forces [48], where numerous references to the literature are given. 

Several theories have been developed which describe the interaction between 
the nucleons by a field of limited range. This field is associated with a 
particle, a meson, the mass of which is between the mass of an electron and 
of a nucleon. Charged particles of intermediate mass have been observed in 
the cosmic rays and lately in rays produced in the BERKELEY cyclotron. The 
simplest assumption is, thus, that the nuclear forces are associated with virtual 
transitions of single, charged mesons causing an interaction term which is linear 
in the meson field function. This gives, however, in the first approximation 
no force between similar particles whereas in fact at least in the 'S-state the 
interaction between two nucleons is almost charge-independent. In the current 
field theories neutral mesons play, therefore, a very important rdle. Hither the 
field is described by neutral mesons alone (neutral theory) or by a combination 
of neutral and charged mesons (symmetrical theory). As to the field theories, 
we would refer the reader to WENTZEL, Quantentheorie der Wellenfelder [43], to 
surveys given by Pavri [41] and Wentzex [47] in the Reviews of Modern 
Physics, to Pauxt, Meson Theory of Nuclear Forces [46], and to ROSENFELD’s 
already quoted book. 
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Since direct experimental evidence in favour of the existence of the neutral 
meson is so far lacking, it may be of interest to try to develop a theory 
of the nuclear forces built only upon charged mesons. Such a possibility is 
given by the meson pair theories possibly combined with a charged linear theory. 
In these theories the nuclear forces are associated with virtual, simultaneous 
transitions of pairs of mesons, one with positive charge and one with negative. 
Thus we get an interaction term which is bilinear in the meson field function. 
Starting from Frrmi’s theory of the f-radioactivity analogous theories were 
developed previous to the meson theories by WerntTzEL [37], Gamow and 
TELLER [37] etc., using pairs of electrons. Then several meson pair theories 
have been presented. Wentzet [41, 42], assuming scalar mesons (spin 0) and 
the simplest spin-independent interaction term between the meson field and 
the nucleons, calculates the interaction potential between the nucleons both 
rigorously and in an approximation of the second order perturbation. JAUCH 
and Loprs [44] have calculated the potential assuming a spin-dependent inter- 
action, scalar mesons and perturbation theory. MarsHax [40], MarsHak and 
WersskopF [41] have in the perturbation approximation made calculations as- 
suming spin 1/2 mesons, while Jaucu [42], Hourret [43], WigNER, CRITCHFIELD 
and TELLER [39], CrrrcHFieLD [41], CrircHFIELD and Lamp [40], have given 
rigorous treatments. KLErn [43], Bruin and Hsatmars [45, 46] have treated 
the case of vector mesons (spin 1) with some different spin-dependent inter- 
actions in the perturbation approximation. Rigorous treatments of spin-inde- 
pendent interaction terms and calculations in the strong coupling approximation 
assuming spin-dependent interactions have been performed by Paui and Hv [45], 
assuming both scalar and vector mesons. In all these cases only charged 
mesons have been used. However, pair theories with both charged and neutral 
mesons have been presented by Noma [46,48] (spin 0 and spin 1/2 mesons) 
and Bruxtn and Hsatmars [47] (spin 1 mesons). As a result of the general 
field theory developed by Kirin [46] some interaction terms, bilinear in the 
meson field function, will occur. 

The purpose of this paper is to give a more complete discussion of the 
charged scalar and vector meson pair theories of the nuclear forces. We shall, 
however, limit ourselves to a treatment of the problem in the perturbation 
approximation. In a paper, published at the same time as this, Hyatmars [49] 
deduces the expression for the contribution to the HAMILTONIAN, corresponding 
to the possible interaction terms between a pair of scalar or vector mesons 
and a nucleon, fulfilling the requirements of LorENTz-invariance and charge- 
independence. In Chapter I of this paper the corresponding nuclear interaction 
potential is calculated in the perturbation of the second order, neglecting the 
recoil and assuming the nucleons to be fixed during the transmission of the 
field. Even the terms depending on the velocity of the nucleons are calculated 
in order to get an estimate of the errors which arise when these terms are 
neglected. The calculation is analogous to that used by Kuern [43]. In Chapter 
II the question of convergence is treated in some particular cases, assuming 
that the nucleons have a finite extension. Further some criterions for the 
validity of the perturbation calculation are deduced. In Chapter III, which 
will be published later on in this journal, a quantitative analysis of the proton- 
proton scattering is made assuming a particular choice of pair interaction and 
of meson mass. The analysis would, however, be about the same with a dif- 
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ferent choice of mass and interaction term. The problem of the deuteron and 
of the neutron-proton scattering is only qualitatively discussed. As distinguished 
from the former case, forces arising from the linear charged meson theories, 
may play an important réle here in addition to the pair forces, and thus, 
these cases are not so decisive for the possibilities of the pair Mheoness 

Finally I wish to express my sincere thanks to my teacher Prof. 0. KLEIN 
not only for having proposed these problems to me but also for his helpful 
guidance and inspiring instruction. Further I wish to thank Fil. lic. S. Hysat- 
MARS for close cooperation on performing this work, Fil. kand Mrs. B. Mosacx, 
who has carried out. most of the numerical work, and Mrs. W. Cameron for 
looking through the English translation. To the SwepisH Atomic ComMITTEE, 
from which I have received grants during the course of this investigation, I 
should like to express my gratitude. 


Chapter I 


Derivation of the nuclear forces in the perturbation of the 
second order 


§ 1. The Hamitronran 


Letting an index s indicate scalar mesons and an index v vector mesons, 
we use a scalar U;(r) and the three space components U,,(r), U»,(r), Uv, (r)* 
of a four dimensional vector as the field functions of the scalar and vector 
meson fields respectively and put P;(r) and P,(r) as their canonical conjugates 
satisfying the ordinary commutation rules: 


[Ue(r), Pe(r')] = [UN0), PX(r)] = éhd (r —r) (1 
. [Ur (r), Pou (r’)] = [U2 (r), Phe(r’)] = 6h by (r —r’) (2) 


all other commutation relations being zero. The fourth component of U», (r) 
can be expressed by the P,,(r):s in the approximation we shall use. 

The wave functions of the nucleons y,(r) (v= N indicates neutrons, » = P 
protons) fulfill the JorpAN-WIGNER conditions of antisymmetric quantisation. - 
We put (summing for y= WN and P) 


r) = av (r) ps (r) 
r) = Pugh (r) oy, (r) 
= PENG r) ip» (r) (i= 1, 2, 3) 
Si(r) = aA Jeop(r) (=1, 2, 3) 


1 Or as a vector U(r). Space vectors are denoted by fat types. 


1 BS) 


OLOF BRULIN, Meson pair theory 


where 64, 02, 03, and 6 (components 04, 02, 03) are the ordinary Dirac ma- | 
trices. 
We put further the meson mass= pm, wc/i =x; the nucleon mass = M, | 
neglecting the mass difference between the proton and the neutron. | 
We write the Hamrttronran of the system of meson field and nucleons 


and 


Hie es (49) | 


where the index 0 refers to the unperturbed meson field, the index n to the © 
nucleons and the index 7 to the interaction between the meson field and the 
nucleons. 


We have 
Hn =— >) | vi (r) (a1 6: pe + e3Mc*) y,(r) dV. (5)i 


As, however, the mass of the nucleon is large compared with that of the 
meson, we shall assume that we may neglect the motion of the nucleons (cf. 
e.g. Kern [43]). Thus we shall disregard the time-dependence of the wave- 
functions of the nucleons. 

Further we have (see e.g. Hyaumars [49]) 


i [set Pal Pan a aa be Us (r) Us(r) + 


+ grad U;(r) grad Cir) aV (6) 


and 


82c : ? ; : 
Hy fi as (x? Py(r)- P35 (r) + div Py(r) div P3(r)) + 


Laz 


1 ‘ 
p alte Us(r)- Uy(r) + rot Us (r)- rot uriryl ay, (7) 


We expand the U:s and P:s in a Fourier-series, whose functions are periodic 
in a cube of volume V, where V tends to infinity. 
In the scalar case we put 


V4 
i (3) 
OPT 2 (kh) (n* (k) — E(k)} et 
Pe(r) = 7S A(K) kn (R) + E(B} et 
* 1 (9) 
PS (r) = Tr Da A(R) (E*() + (Be 
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and 
i ; 
grad Us (r) = Dy ke (h) {E*(k) — n(h)} (10) 
where 
u(k) =V4ach (2 + x)-* (11) 
rm) = |/ a+ a4 12 
a l6ac. em} eo 
and the following commutation rules are valid 
[&(k), &* (k’)] = [4 (k), 1* (h')] = Onn’ (13) 


all other commutation relations being zero. Introducing these expressions into 
(6) we get 


Hos = ch > VRB + 2? (Nk) + Nb +1) (14) 
where 
N* (k) = &* (k) E(k); (15) 
N~(k) = *(k) 9 (k); (16) 
N*(k) and N (k) are diagonal matrices with 
(N° | Ne N ) = N° = 0 (17) 
(N |N-b&|N)=N 20 (18) 


where NV* and N are integers. The only terms of é(k) and é*(k) or of 7(k) 
and 7*(k) different from zero are 


(N*|é(|N*+1)=VNT+1 (N+ 1/e*|N*)=VNTF1 (19) 
and 


(N|y@|N + =VN+1 (N+ 1|y*®|N)=VN +1 (20) 


where we have disregarded the time dependence. Our unperturbed system can 
thus be described by N’(k), N (k) and some quantum numbers, representing 
the nuclear system, which we do not need to specify any further. As we are 
interested in the lowest quantum state we shall put both 


iN erand. a\' == ()) (21) 
In the vector case we put 
Us (r) = —= >, Dd) po (k) &» (k) (EX) — ny (de) ef 
VV k v=1 
(22) 
Ut (r) = = > Dd wilb)a(h) (gh — Eb) 
V hp 7yc—ik 
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Pat) = os yy Ay (Kk) &» (Ke) (no (he) + &,(h)) e” 


VV ko v= (23) | 
PS (r) =D, > Aa le) es (HY (G20) + (he | 
V ko v= 
1 * \ thy 
div Up(r) = ve >) pe (Hb) Bee & () (ERB) — p(B) e* (24) 
k v=1 
rot U(r) = => fbr (ke) kk X &y (ke) (E5 () — (he) (25) 
k v=1 
Sheek Wales Vi | > Aa (Ke) 85 (le) (3) + Ey(do) eo” (26) 
k v=1 


where! ¢,(k) are three mutually perpendicular unit vectors of which we chose 
&3(k) parallel to k i.e. 


k = k &3(k) (27) 
an ae 
ia (k) = fa (k) = as sad ak (28) 
3 (k) = Vienne ace 
and Es 
=e = SE de 
is (8) == ta(h) = |/ gin + 08 
oe (29) 
Sal Joa! ayy 2 ee 
A(h) = | 6x i (Ra He) 
The following commutation rules are valid 
[&. (he), &" (e’)] = [yy (Ke), 90" (he) ] = Sv" Orca (30) 
all other commutation relations being equal to zero. 
Thus we get 
Ho =ch > SVP +2 (Ni) + Nv) 41) (31) 
kk v=1 
where 
Ny (k) = & (k) & (k 
. : (32) 
Ny (k) = 15 (k) 7» (k) 


The matrices N;(k), N,(k), &(k), &(k), 7, (k), 73 (k) have the same properties 
as the corresponding matrices in the scalar case (17)—(20) and we get the 


corresponding quantum numbers. As before we shall only consider the case 
NSS N 5 ==0: 


* A point denotes a scalar product and a cross a vector product. 
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If we consider all interaction terms between a pair of scalar mesons and a 
nucleon, satisfying the requirements of LoRENTz-invariance and of charge-in- 
dependence we get (cf. Hsatmars [49}) 


2 
He = aa i) dV - (33) 
Cs S3(r) Us (r) Us(r) + (S;") 


82¢ 


+ Cx S3(r) ie grad U*(r) grad U(r) —( : ) Pye r) + (Sy) 


C 


Prarie 
eC. [2 So(r) (P(A) Ur) — U* (ry) P¥() + 


a Siric 


M4 


(grad U*(r)- U(r) — U* (r) grad ur) =f 
m 


8 a 
Os |= So (r) (Pr) U(r) + U*(r) P¥(r)) + 


+ S; (r)— (grad U*(r)- U(r) + U* (r) grad Ur] + 


+ iCss [sim 2 grad U*(r)X grad U(r) + 
| (S3’) 
Fa ag aes (Fd ee otad.U (7) —.eradgU Fin) F* |} 


where the C:s are arbitrary real constants, and for the moment the index s 
is neglected in Us(r) and Ps;(r). This expression, however, is only an ap- 
proximation, correct up to and including the first power of C. We put 


(&° + x°)~* = (hk) (34) 
ae (35) 
an (", Rl) == 9" (H") BW) + E(R") 9 (Rt) — a ("9 (le) — & (BY BB) 
(KB) = gh (SR) + ECR) 0 (R) + RO) n() + ERE) 
B (kK, k’) = * (hk) E(k’) — E(k") 0 (h’) — 78 (B") 9 (B') + E(B) oF (B’) 
B (hk, ke’) = 1 (hi) E(k’) — E(k) (h') + * (B"') 9 (BR) — 8 (B") &* (B') 
Introducing (8)—(12), (84)—(36) into (33) we get 
4 ; yn Cel Oe 
‘i Os Sg (r’) w(B) (B') a (k", ke’) — (Sy') 


{ 
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— 0% 83(r’) |  B(kK) uw (k') (ki, R’) + A(R) ACh) a (B", | - 

+ Ce [s. (r’) (AK) wk’) Bik", kb) — wk Ak’) BR", k’)) — 
—Sy(r') 2 ("+ Ry (hy (hack, WY] 


4 Oe |- So(r’) (Ak) wk Bk", kh’) + wlh'AK') Bk’, k)) + 


Sir) (YB) (lB) ("| 
yO E Se (r’) p(B” XR) ju (Ke) je() a(R", k’) + 


+ S,(r’) : (hk! ACh) a(k’) Bk", kh’) + ke whe) Ak) BR, &)|I. 


Similarly we get in the vector case 


JU C 


{Ci S3(r) (us (r) U(r) — a (222) div P(r) div P* r)] ~ 


) Pw BAe) = 5 rot U*(r) rot U(r)] a 


8 7¢ 


+ 4Cre [s. (r) 3 (Pr) U(r) — U* (r) P¥ (nr) — 
te (r) 2 (rot U* (r) x U(r) + U*(r) Xrot U(r) + 


82¢ 


2 
ae ("2") (Pir) div P*(r) — div P(r): P* »))| As 


+ Cre [S007 —— (P(r) U@) + UFG@) P*¥r) — 


194 (5) : (rot U* (r) X U(r) — U* (r) X rot U(r) + 


8 2 
- ( a (P(r) div P*(r) + div P(r) P* r)| 


2 
an oe 8 2 I 
+ 71Org [ ssi ((°2") P(r) X P* (r) — 2 rot U*(r) X rot U(r)] = 
LCS Cpe : : = 
== 95\(77) pe (P(r) Xrot U(r) + rot Utinx Ptr) | + 
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(V3") 


1 : 

mer (rot U*(r) div P*(r) + div P(r) rot Uir)] ae 
1 

+ 8, (r) 7 (rot U*(r)- U(r) + U*(r) rot uir)| — 


8 2c 


+ i0v | s (r) (Pr XU) + UF) XP*(r)) — 


1 tone : 
= 2 (rot U*(r) div P*(r)— div P(r) rot U'r) 


iT) d 


a 


) F (20) 


rot U*(r)- U(r) — U*(r) rot ur)| as 


; 8 ; 
+ Cs So (r) 5 (P(r) rot Utr) + rot Ue) Per} (V5) 
uw 
where the C:s are arbitrary real constants and the index v in U,(r) and P,(r) 
is neglected. As above, this expression is correct up to the first power of C. 
We put 


(39) 


Bare: (He, ke) = nh (RY) EF (Ie!) — Ey (le) ee (Re) (He) ge (Be!) — Sv (Be) (Be). 


Introducing (22)—(29), (39), (40) into (38) we get 


ye Me Kk on” 


ror =1 
| 850!) (— pr (Be fe (8 (60 6) er (RR) — 


1 a _ thes / vt / £ = vt , 
Aa 8") Ay (HE) («00 (Be = 6 (H9) rer (BB ) at 
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=z Ch S3(r’) (i (k’’) Ay! (k’) (€- (ke) : &;' (k’)) Cy" y" (Kk, k’) oe 


5 Shir (K’’) prt (he!) (he! X & er (B) (Be X € )) er (B, K)) + 


+ Ove Ex (2. (He’") fay’ (B’) (&r'r (He) > & Ue’) Breve (Be, ) — 
— fy (Ke) Aye (Ie) (8 (Be!) > pe) Brrr (He, K)) “ (Vo) | 
At +(e”) ft (le) { (Bo X 84 (B') X &x (’) — 

— & (BX (BX eh) 5 ape (RB) + 


4+ Aye (Bel) dye (Ie) Lepr (Bel) (Bel 6p (He) (Ie > 8 (Te!) 8 (BY) ett (Be, K))| es 
10%, |S (r’) (— J (He”) fart (B’) (Beer UY 89°(6) Bret (Be, Be’) — 


— pip (Kel") Aye (Be!) (8p (Bel) > & (Be) Brrr (he, K) + 
if (V2") 
Sar) (ee py" (Be) pay’ (Be) (ke X 8p () X &y (B) + 


ar Ey" (k’’) »* (k’ XE," ‘k’))} Ay!" y! (hee k’) =a 


ir: (Kh) Ay (he) ep (Be) (Bl & he) — (e+ & B) & (B’)} Gyre (k, K)| a 


is iCal Solr a (k') Ay (K’) (pr Be") X &y (he’)) Gyre (ke, he’) + 
4 fy (RY) (Re) <a {(R”X ey) X (BE X 6B) (K",#)) 53 

. (V3’) 
. ae Ar (K") far (Ke!) Ls (Be) X (Bl Xp (!))} Brot (he, Be) + 
+ Mr (k’’) Ay (k’) (kX ep Bh) X €r (B’)} Brrr (RY, m) | + 
+ 40d3 | s; (r’) (- Mart (B") tte (h') (8 DX €p 2) ape (BM, K)) te 

(Vs) 
- Sar) (Ben (h) Hr (h’) (C= &y(")) &e (B')} Byrn (B", Be’) + 

Dy (k! 


tH ptats (el) Ay (Be) {50 (0) (> 890(B)) | Brrr (Ie, om) if 


EO, | s. (r’) (— Ayer (He!”) jure (Ie) ("8 X (8) Byer (Bl, Bel) 


les ~ > 
a pe ee” (k ) Ay’ (k’) (kh X Ep! k'")) (k’ + Ey (k’)) Br''r! (k”, k’) == 
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tor (Be) Are () (E(B) X 0 (")) Brrr (Be, Be’) + 
ibs . (Va) 
Fwy A (He) te (He!) (e+ 6 (le) (eX eye) Brey Re", BD) + 
ilar 
+ 8, (r’) a (i (he!) tr (Be) (he! X &p (') + & (’) — 
— (be) + (kX ee (B))} ane (k", m)) | + 
+ One | s. (r’) (2. (K’’) tr (h’) (80 (kX & e)) Brie (ke, ke’) — 
ae = md 
a beri (He) Aye (Ke) (eX & pe) (Bel > & er e)) Bore (", k’) + 
tee (Ue) dee (Re) (6B) X 6-1 HY) Bor (Be, Be’) (V4) 
1 he , wr 
Saber (Be (RO) (Ry Re) (BX £7) Bre (",B')) 
i / is = ite eS , ft 
P+ (7) (= ae (he’") pert = (he!) (BX & (Be) &e (he) + 
1 + ey (ke) (kX en (k))} ape (bk, m))| + 
, if = = 
+ Cis S:(r’) = (i (K") jay (I!) 0 (le”) - (lt X °C) Brrr (i, Be!) — 
“ ~ ua (Vs’) 
bet (BE) Ree (Be!) (BEX Er Be!) + &e (BY Brrr (B", K’))]- 
) § 2. The perturbation calculation 
A. Vector mesons 
We will begin with the vector meson interaction as the arguments, applied 
to this case, can easily be simplified so as to be valid in the scalar case too. 


As the perturbation of the first order does not give any interaction energy! we 
may limit ourselves to the perturbation of the second order. The diagonal ele- 
ments of the matrix of the energy of the perturbation of the second order are 
sy (| Hi | 0") (n" | Hin’) (42) 


CARAFE aa Ent — En" 


where the quantum numbers are given immediately after (32) (cf. (21)). 
The matrix H; is given in (41) (or (37)). As in the quantum numbers n’ in 
(41) (or (37)) N+’ and N~' are put =0 and thus N+” and N~" are = 1 (see 
(19) and (20)) the only terms in the matrix product (42) which are different 
from zero contain a factor (cf. (32) or (15) and (16)) 


En (kh) ne (Ke) ns (U") Es" (U) = Ov 5 Orr's’ One Ont’. (44) 
1 The neglected C? terms in (38) (or (33)) give in some cases contact interaction terms. 


Such terms may be neglected, however, owing to the break down of the theory at small 
distances (cf. Chapters II and IIT). 
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From (31) (or (14)) we have 
Bit = Env SS h Cc (Vie? + x“ =e Vk? + x2). (45) 


From (44) and (40) (resp. (36)) it follows that in (42) we may put a= @ and 
B= and we get 


cs’? 5 (UU) yr (he, Be’) = Ott" Ost” Our! Ov” (46 a) 
ats’ s' (LU) Byte (BY, i!) = Og v” Ost r”? Ov’ Orr” (46 b) | 
Botts (UU) ayrry (Be, Be’) = — Ogre Os’ Oe! Orn" (46 c) 
Bots (U,V) Brn (Ke, B) = — Sarre Ostet Or! Or'e". (46 d) 


| 
| 


The interaction energy AH is obtained if in H; ((41) or (37)) instead of the © 


wt 


indices k”’, k’, 7’, r’ we write Il’, I’, s’, s’ and multiply this expression with 
the original H; where, however, instead of the variable r’ we write r”’. Ac- 
cording to (46) I’, I’, s’’, s’ are changed into k’, k’’, s’, s’’, and in the case 
of (46 c) and (d) even the sign is changed. In practice this means that to get 
AH we have to multiply H; (41 or 37) by itself, having in the first factor 
exchanged k’’ and k’ and 7’ and 7’, in the second factor changed the variable 
r’ into r”’, and replaced the products of a and f with a plus or minus sign 
in accordance with (46 a,b) or (46c¢,d). As in the sums and integrals the 
notation of the indices and variables is immaterial, we see that, if we change 
the notation suitably, many terms in AH occur twice. These terms have 
different signs if one term has af as a factor while the other has Pa (46 
b and c). Thus a great many pairs of terms cancel each other and other terms 
double themselves. It is, particularly in the cross terms between different inter- 
actions, often useful to exchange the indices in some terms and alter the signs 
of k’ and k”’, «&(k’) and &,(k’’). This means that the sign of the term is 
altered, if it has an odd number of factors k. In this way, we shall again 
have some terms which either cancel one another or double themselves. 

We shall now let the volume V tend to infinity replacing the summations 
with respect to k’ and k” by integrals. In accordance with an earlier paper 
by Hsaumars and the present author [44] we shall multiply the integrand by 
a function f(|k’|,|k’’|,a) so that the integrals converge. a is a parameter 
(cf. also Noma [48]). In order to evaluate the k-integrals we will introduce 
polar coordinates taking r’’—r’ as the polar axis. We shall use a rectangular 
coordinate system too, the z-axis being the same as the polar axis. The com- 
ponents of k in the polar coordinate system being (k, 3, ~) we may put, 
using (27) 


&, (k) =z sin p — y cos 
& (k) = % cos 8 cos y + ¥ cos Y sin m — Zsin 0 (48) 
£5 (k) =a sin 9 cos + § sin D sin @ + 3.¢08 9 > 


where z, y and 2 are the coordinate unit vectors of the rectangular system. 
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We shall put 
|k'|=xp, |k’|=xq, cos # =u, cos 0” = v, |r’ —r'l=r (49) 


and shall use rather often the following elementary integrals 


20 27 27 

[cos pdp = | sin p dy = [ cos p sin pdg = 0 (50 a) 
0 0 0 

27 27 

| cos pdy = | sin? pdyg =a (50 b) 
0 0 

27% 

[dp=2n (50 c) 
0 


From the integrations with respect to gy’ and ” and from the arguments on 
the preceding page it follows that all terms of the interaction energy of the 
second order perturbation arising from different terms in (41) (or 37) are zero 
with exception of those arising from V, and V,4i or S, and S,4+1. 

We get from (41) 


1 e 2 ar 
AH. (a) =—5'3() we { favav’ | [avag: 
0 0 


St mt pt 
as fo, qa)r ¢ is / | OND A ad (51) 
Vp? = 1Ve? + 1(Ve? +14 Ve+ WJ} r, 
- {On Ss(r’) Ss (r’") (6 + 2 (p? + g®) + 3p? Q? — pp? (vu? + 07) + BpPGurr® + 
+4pquoVp? +1Ve4+1) + (Vi1') 
+ On S3(r’) Ss (r'’) (6 + 4 (p? + 2) + 6p? g? — 2p? G2 (u2 + v2) + Op? gurv? + 
+ 8pquoVp? +1Ve4+1)+ (View) 
+ On Ox Sg (r’) g(r”) (— 12 pquv — 12 Vp? + 1Vq@? + 1) + 
+ Cre [Sp (r’) So (r'”) (12 + 6 (p? + a”) + 2p? gq? + 2p"? (u? + v*) 6p gu? vo? — 
—12Vp? +1Vg? +1) + (Vin") 
mt 9, (r’) S;(r’) (5 (7 = 0) (p? uw? + gv) +477 gt 4p? Gurv® + 
+ 8pquvVp? +1V_e +1) + (V2'2') 


‘ Deo Ay Oa ae Oe 
+ Siz (r’) Siz(r’’) (— (p? + a?) +3 (p? v2 +g? 0?) + 2p? Gg? — 2H? ¢? (u* + 0°) 


— 6 prgurv? —8pquoV p? + 1V@+1+ 
Si 12 pquv)] sl 
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+ Osa [So (r’) So (r’”) (12 + 6 (p? + g?) + 2p? q? + 2p? Gg? (u? + v*) —6 pg? ur v® + 


+12Vp? +1Vqg? +1) + 


+ §,(r’) S (r’’) (5(p? + @)—(pwt+ gv) t+ 4p? @t 4p Purr? + 


+ Spquv Vp? + 1Vq? +1) + (Vera) | 


+ Siz(r’) Siz (r'’) (— (p? +.@2) + 3 (p? u? + g? vo?) — 2p? @? (uv? + v?) — 6p? Pwr vr— | 


—8pquoVp? +1Ve+1— 
= 12 pquv)] aie 
Gls LS (r’) S3(r’’) (4 + 3(p? + @) — (pw? + ov?) + 4p? PF — 4p gu? v?) + 


+ Ssz(r’) Ssz (r’’) (— (p? + @?) + 3 (p? u? + g? v?) — 2 pq? + 2p? Q? (u? + v7) + 


+ 6p? g?u2v? + 8pquvVp? + 1V_q? +1) + (V3'3’) 
+ Sp(r’) Sy (r”) (3 (p? + ?) — (pu? + g?v*) + 4p? g? — 4p? q? uw? 0?) + 
+ See (r’) Soe (r”) (— (p? + g?) + 3(p? ue? + G?v*) — 2 pg? + 2p? g? (u? + v%) + 
+ 6p? grurv? + 


+ 8pquoVp? +1Vq? + 1)) + 


+ Ova [Sg (r’) Sa (r’") (4+ (p? + @?) + (p? w+ Gv?) +p? @? (u? + 0?) — 2p? Pu? v®) + 


+ Ssz (1) Sse (r’’) ((p? + @?)—3 (p? u? + gv?) + p? ?@—2.p? G2 (u2 + v2) +3 p2gru2 v2) + 
+ Sz (r’) Sp (r”) (2 (p? + ?) + 2p? q? — pg? (u? + v*)) + (Vas) 
+ Sez (r’) Sz (r’’) ae 2p 2g + 3 pq? (u2 + v7) + 4pquoVp? + 1V¢q@? + 1)) + 

+ Cvs Crs [S53 (r’) Sy (r’) (— 8 Vp? + 1 Vg? + 1) + 

+ Sz (r’) Ssz (r’’) (— 8pquv) + (V3'3"") 
+ (S3(r’) Sz (r’”) — Sez (r’) Soz(r’”)) (8 pquv)] + 

a3 Cra [So (r’) Sp (r”) (8 +4(pP? + @)+ 4p P+ 4p? Q@uiv® + 


+ SpquoVyp? +1Ve—i +1— 


— 8pquv—8Vp? + 1V_e +1) + 


+ Soz(r’) Soz (r’’) (2 p? @? — 2p? @? (w+) — 6p gue? — 


—8pquoVp? +1Vg +1 4 (V4) 
+ 8pquv) + 
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+ 8, (r’) Sy (r’’) (4 (p? + P) + 2p + Qp%Q? (u? + vw?) — 6 pr guro? + 
35 Spquv)| Te 


+ Ors [So (r’) So(r”) (8 +4 + 2) +42 @t4pewr + 


+ 8pquoVp? +1Ve+1+ 


+ 8pquv + 8Vp? +1V_e2 +1) + 


+ Soz(r’) Soz (r’’) (2 p? g? — 2 p? 2 (u? + v2) — 6p? @utv? — 


—8pquoVp?+1Ve_+1— 
— 8pqur) 3F 

+ Sy (r’) Sy (r'’) (4 (p? + @?) + 2p%Q? + 2p? Q? (wu? + 0?) — 6p? gtu2o® — 
—8pquv)] + 


O45 So (r’) Sa (r’’) (4 (p? + g?) + 6p? g? — 2p? @? (u2 + 0?) + 6p gure? + 


+ 8pquoVp? + 1Vq? + 1) + 
+ 4 Om Ove (Ss (r’) Siz (r’’) — Siz (r’) Sg (r’)) (4 (pu —qv) + 


+ 2pquv(pu—qv) + 2(pu—quv) Vp? + 1Vq? +1) + 
FF On Cre (So (r’) Saz(r’’) — Siz (r’) Sq (r’”))-(— 2 (pu — 0) — 


—4pquv(pu—qv) —4(pu—qr) Vp? + 1Vq@ +1) + 
+ Ove on [(S, (r’) Sez (r’’) — Sez (r’) So (r’’)) (ae. (pu —qv) + 


(Vara) 


(Vs's') 


(Vi2") 


(View) 


+ 2p? qu—2pgut 2pquv(pu—qv) + 4 (pu — qv) Vp? + 1Vq@ +1) + 


+ (S,(r’) X S3(r’) + S3(r’)X S,(r’))2-(— 2 (pu — qv) — 


—pqvut peut 3pquv(pu—qv) + 2 (pu — qv) Vp? ab 1V¢ ae 1)] ae 


+ 4 Ore Cvs [(Sq (r’) Sze (r”) — Sax (r’) So (r’”))- (— 2 (pu — av) + 
+ 2pquv(pu—qv) + 2(pu—qr) Vy? 1 Veeck Le 


+ (S,(r’)X S3(r”) + S3(r’) X Si (r'"))e: (— 2 (pu — gv) + 


+ 2Ipquv(pu—qv) + 2(pu— qv) Vp? + 1V¢q@?+1))+ 
Boe Ons Cos [(S3z (7) S, (r’”) — Sy (r’) Sez (r’’))(— 2p? qv t+ 2pG@u— 
—2pquv(pu— qv) —4(pu—qu) Vp? =. 1V_q@ ae 


+ (S2(r')X So (r”) + So(r’)X Sa (r”))e-(— 2 (pu — av) + 


+ p2qu— pgu—3pquv(pu—qv) —2(pu— qv) Vp? + iV¢@ + 1)\+ 


(Vo'3’) 


(Vers) 


(Vs) 
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+ § Cog Coa [(S32 (r') Sy (2) = Sy (Se (P')) (+ 4 — 90) — 
—I2prqut+2pQ@ut 2pquv(pu—qv)) + (Vara) | 
+ (S,(r’)X So (r’’) + So(r’) X Sa (r’’))2- (2 (pu — qv) + 

+ pqvu—peut pquv(pu—aqv) + 2(pu— qv) Vp? + ie Ll 
+, Osa Cus (Soe (r") Sa (r") = So (Sent) (ipa — one | 
+ 4pquv(pu—qv) + 4(pu—aqu)Vp? + 1V_@ + 1} (Vis) | 


all other interaction terms being zero. 

In this equation we put xr=2 in xrpu and xr=y in xrqv. The func- 
tion /(p,q,a) in (51) is assumed to tend to zero as p or q tend to infinity if 
a-+0, so that the integrals in (51) converge uniformely with regard to x and 
y through the entire space. ot: 

Further 


ae f(p,@, 4) = 1. (52) 


We might e.g. put / (p,q, a) = e% Pt, 

If we write the integrand of (51) as derivatives with regard to x and y it 
is permissible to remove the differential operators outside the integral signs as 
the integrals are uniformely convergent. 

We shall use the following integrals 


Li(@y.a)= 
ae pe dpdat(p,qa rae 
=tev| = Bae dudv ei @pu-yn) — 
Vp +1VG+1(Vp? ST. +1), 
0 0 
a= al ——— EE ERS = sin'2'p sin (53) 
J) Vor + 1Ve + 1 (Vp? +1 + Ve? + 1) Z° 
and 
I, (x, y, a) = 
eda dare (p,q, a 
=a | | dudv et @pu- yqv) — 
Vo? +1 4+ Ve + lJ 
af pqdpdgqft (p,q, a) dene a 
Vp +1+Ve4+1 a 


Then we can write (51) as derivatives of I, (a, y, a), Io (x, y, a), eit (x, y, a) 


1 
de las i 
an a o(@, Y, a) 1. ¢. 
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~ +1 +441 
an gq’ dp dat (p,q, a) J ae YQ). 
; Vp? se eater: Ta ott) 
ye 
ae hae 
i qgdpdat(p,a, a) e| Piibena 
Vp? ae Oe ee Va eat) : 
Sib aul 
pt i : pqdpdqt (p,q, a) ye cin a 
LY Vp? +1Vq@? + 1 (Vp? +1+V@ +1) 
0 
We GE 
gy age (ZY; a) (56) 


and correspondingly for the other terms of (51). 
Thus we get 4H,(a) expressed as derivatives of (53) and (54) 


AH, (a) = ~3(<) uc ff avay”. (57) 


ie Se (r’) Sar” {2 ‘ ee | ae eo a 1 @ ets 


Aes y\0xe dy? cy Oa dy? Ox? xy Oy" 
Ce: ieee aS) Cae wal 
i Ay? Ly 73) eae cy EK (Z, Oe G)i-P et 
Ne , ” f 6 0° | 0° 6 os se 2( 0 1 o | yy 
B Coy Bs (r') Sa (r (j2-s cy\0x22 | OY? * cydx2dy “\daeaydy te) 
eo uc \ 5, OF 3) -1| OF an: 
' OY xy +i) 0x72 0Y" a (2, y, a) any LY ICRP Fa 
”  r) Se Or” Caley Di + (Vin) 
| C,, C,, 8s (r’) Ss (r’”) Mon eo) 2 (2, Y, a) (Vin 
he ; ete Oe EO ee Os \y piece © Oe (= Leman 
me Exc ee (fe alas ml ' cy 027 0y? O22 zy OY - 
a ot 1 
| : eo, , ae 
Oy" xy aa 0x22 0Y" =} 1 (@, Y, a) 
— = te(e,%0)) + (V2) 
“LY 
; Re (eis 0 -B00e | se On Vel 4 4 ot 
+ Sy (r’) Sy (r (\- (ae xa) (5 Oyr)ay «cy O20Y 
ERAN ie 
eee NE ir oye) | oe 
EE 7 (enh a) Ox Oy vy (x, y, a) 
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We shall define AH as lim AH (a) having put the differential operators 
a=0 
outside the integral sign, as in (57), so that the integrals left are convergent 
even when a=0. We get from (53) or (54) 


I, (x, y, 9) = lm J, (a, y, a) = 
a=0 


pdpdg sin xp cos yq 


0 
= hm | ae aa P,d, a 58 
a=0 an Vp? t+ Ve +1(Vp? +1 + Ve? + 1) oe 
where 
fH — pdpdg sin xp cos yq (59) 
J V; meen > + 1 eee) 
) 
= [' cos yqgdq ie sin a ee lig sin xpdp ee yqdq 
Vee yey oe Vp? +1 ieee 
0 0 


(taking the principal values of the integrals) 
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(taking the principal values of the integrals) 
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(x,y) i PVE LI 
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0 F (x, y) “cos (w + y)@ 
ESS aK ety). 
0a" | Vg ay 


Using the well-known differential formulas for the HanKe functions Ko (2) 
and K;, (2) 
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0 0 
By et RO ee (63) 
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7 ida f Gea tg 120 : 1 
+ Co Coz (S3(r’) Sy (rn) — Sy(r') 837") 2 |- (S a =) K, (2) — 


+ Cy Cos (So (rr). Sor’) — Solr) So@ 9 4 SiO) Sar Oe ie Sy (r'’))z° 


ess a) Ea) + Ke@l+ 


+ Coe 0% (So(r’) Sy (r") — Sg(r’)Sy(r”) + Sy(r')X Sgr) + ae (r’))e- 


[2 $)e0 Sao] 0 


as Ce Ona (S3(r’) Sy rs) > S) (r’) S3(r’’) ae Ss (r’) x So Yr’ oy = Soir) So r’))e 4 


[epsom] 


+ O%3 Cvs (S3(r’) 8, (r'’”) — Sy (r')S3(r') + So(r')K Sor”) + Sor’) X Sa(r)) 


(Vs'4") 
40 10 20 
[(2¥) xe 8)aa] + 0-0 
+ Org Crs (So (r") (Se (r’”) — Se (r') Solr”) z lle i | Ky (z) + “i Ko |} 


(Va5') 


B. Scalar Mesons 


The calculation is in this case analogous to the vector case. The equations 
(46) are still valid, but the indices r’1r’’s’s’’ disappear. The same coordinate 
systems are introduced where k is given from the third equation in (48). 
Thus we get from (37) 


AH; (a) =— 5 1a(<) me *f favav’ | { avaa (68) 
0 0 0 6 ; 
t ss +1 +1 

PP, 4; pq : ixr (pu—qr) ; 

e dv: 

Wary | | oe 

ayant 

{0% S3(r’) Sy (r’”) 2 + (Srv) 


+ Cs Sg (r’) Sg(r") (2 + 2 (p? + 9?) + 3p? q? — p2g? (u? + 0?) + 


+ 3p? @urv t+ 4pquoVp? + 1Vq2 +1) 4 (S11) 


+ Cs Cs 83 (r’) Ss (r’) (4pquv + 4Vp2 + 1V¢q? +1) + (S11) 
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+ Ce [Sp (r’) So (r”) (4 + 2 (p? +g?) —4Vn®? F1V2 +1) + 

+ Sy (r’) Sy (r'") (p? + @? — (pe? + @v®)) + (Sor2") 
+ S12(r')S1z(r'") (— (p? + @?) + 3 (p?u® + @v®) + 4 pquv)] + 
+ C49 [Sp (r’) Sp (r’”) (4 + 2(p? + g2) + 4Vp? + 1Vqe2 + 1) + 


+ S,(r') Sy (r”) (p? + o? — (p? u? + g? v®)) + (Sor2”) 
+ Siz(r’) Siz (r'’) (— (p? 4 g*) + 3(p2u? + gv?) —4nquv)] + 


Fe ( [S3(r’) S3(r'’) (p? @? (u? + 0?) — 2 p? gq? u? v?) + 


+ S3z(r’) Sgz(r’’) (p? g? — 2 pg? (u2 + v7) + 3 pq? u?v?) + 
f Sa(r) S5(r ) (p? + Q@—(pPwt+ gv) t+ 2p? —pe (w+ v)) = (S3'3’) 


+ Soz(r’) Soz(r’’) (— (p? + @®) + 3 (p? uv? + go?) — 2p? @? + 


+ 3p%q? (u2 + v) + 4pquoVp? + 1V@? + 1)) + 
+ 4 Os C% (Ss (r’) Siz (r’”) — Siz (r’) Sg (r’’)) (2 (p u — qv)) + (Sy So”) 
+ 404 Ove (Ss (r’) Siz (r’’) — Siz(r’) Sg (r’’)) (2 pq uv) (pu —aqv) + (Sy Sov) 


+ 2(pu—qv)Vp? + 1Ve2 +1) + 
+ 1 Ose Crs [(Sp (r’) Saz(r’”) — Sze (r’) Sq (r’’)) (2 (pu — av) — 
—2(p?qu — pau) —2(pu—quv) Vp? + 1Vq? + 1) + 
+ (S,(r’)X S5(r”) + Ss(r’) X Sy(r"))e(p qv — (S23") 
—pgeu—paqur(pu—qr))]f. 


Similarly to (57) we get 


2 é 2 £. vr 
AH; (a) Se len Mere dV’ av”: (69) 
{ "2 , ” 2 g 
3 \Ca S3 (r )S3 (r ee (Ze Y, a) == (S11) 


(2 se Le a4 


+ C,,S3(r') S3(r’”) ie a "py Ox OY? 


Ox? Oy 


Cpa s0 ee | “ | Gat!) Pty cone iS) 
laa ay Oy Oyray Ox?) | ~ Ox? dy vy | ie 
Oo ae 
ote 2 
ai yy heya] 
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(5, pent laie, UN, Orel , 
Ox? LY Oy? dy? Ly 0x2 i 1 (2, Y, a) ie cP ieee «) oi 
+ Os1 O52 (S3(r') Siz(r”) — Siz(r’) Ss (7) (2(5, ae -) Ly; (x, y,a)) + (Sy2") 
\ / Ox dy ry 1 9D) 1/2 
” ’ peg (ok 
+ 1 Cra (850) Stel") — Sta(r)S = ber - (Sv'2") 


ier jx 0 a) al 
a qi 1 
Ox Oy) ay 


pei ue leeko ty 302 \) pao al 
ay LY da? Oe cy OY" 3) (2, y, a) Ae A Voes (2.1 @)) + (S2’3’) 


+ Cs2 Css oe Sez(r’’) — Sez(r’) So(r’’)) ( ( 


/ tae ay ef Onan LOe 
Poy OCS eG). Sar’) XS (rr), eee as 


Benet 7 + Nernst ah 
Oa xy OY? Oa? Oy Oxy?) xy teu) 


where x=y=x|r’—r’| and all other interaction terms are zero. Putting 
a=0 we get (cf. (67)) 


ee se) oye [{favare. (70) 
Aci S5(r') Sgr") (— 2 Kt )) + (Sy) 
+ OSS 0") (— (Br + 4 Pl — (EP + s)he) + Ger) 
+ O91 Cn 83 (r’) 83 (r’’) (- (3 + | K, (2) —3K, ) + (S11) 
“Nei [s (r’) So (r”) = K, (Qe 3K, ©) 2 

+ Slr) Sit") (— 4 Bie) — 3 Kol)) + (Se) 
+ Sie(r') Stele") (2 Ka e) + 4 Kole)) | + 

+ Og] Sie) Su) (— 2 KL — Ko) + 

+ Siz(r’) Siz (r”) (3 + 3 Ki) + = Ko )| i (See2) 


145 


OLOF BRULIN, Meson pair theory 


+ 3 | Scr Sst" (AR + Le) wa) + (B + 3) ol) + 


x 


14059 20 70 2 
+ Sbe(r') Bie(r") (— (He + ha (+ 3) Kole) + (Ses) | 
20 2 10 
+ S53 (r’) So (r”) (= at ) IG @y Z Ko (9) = 
140 20 oe & | 
+ Sez (r’) Sez (r’’) (~ ~6 4 Ky (z) a" (3 1. 5) K )| an 
, tt , 3 1 
oe On Gr (Sar) Sy) — 8; r)S. ir Ne (30 sir +Ko()] aie (S1'2") 
60 9 
ate Ce Os (S3(r’) S; r’’) ae S; r S3(r’’))z (( 5] Kk, (z) i 


i 
at (= a 4) Ko ) oF (S12) 
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Chapter II 
On the question of convergence in the meson pair theory 


§ 1. A rigorous determination of the total field energy in some particular cases 


As is well known, we get an infinite contribution to the meson field, if, 
assuming pair interaction between the source and the field, we introduce a 
point source in the meson field. The interaction energy between two point 
sources is infinite too. Consequently, it is necessary to make another assump- 
tion as to the shape of the source. We might assume a source of finite ex- 
tension. Such a source could for example give a rough picture of a nucleon 
that is, as is suggested by Ktiern [48], assumed to be built up from several 
smaller particles which are the real sources of the field. Our assumption is 
not, of course, relativistically invariant and we shall only use it when the 
source is at rest.’ Chapter I can immediately be applied to a source function 
of this kind. Assuming this source, we want to investigate the self-energy of 
the source and the interaction energy between two sources. We shall at first 
use a rigorous method restricting ourselves to a few particular cases. 


* This is not the assumption of an extended source usually made in the pair theory. 
(See e.g. WenrzEx [41, 42], Pautr and Hu [45)). 
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In the lowest state (the state of »no mesons») the quantization of the field 
gives a total field energy equal to half the sum of the moduli of all the 
energy values which are obtained from the energy states of the field equation. 
As is shown by Hysarmars [49], this is true not only if we have an unperturbed 
meson field but also if we have an interaction with a source which is bilinear 
in the meson field (pair interaction). In order to obtain these energy values 
we shall imagine our system of field and source included in some space of 
finite but large dimensions, requiring our field function to be zero on the sur- 
face of the space. Thus we get a discontinuous energy spectrum. We sum the 
moduli of the energy values letting the dimensions of the space become in- 
finitely large. In doing this we shall bear in mind which terms of the energy 
belong to the unperturbed system, which terms give the self-energy of the 
single sources and which terms give the interaction energy between the sources. 

We shall only consider scalar mesons and we shall assume that the contri- 
bution to the Hamittronran from the interaction between the source and the 
meson field is 


4H=2 | ueu gay (1) 
82 


where U is the meson field, and y2 is the source function." The integrals 
must be taken over the space, defined above. Of course, [aQav ought to be 


finite. We shall, however, discuss some cases where this condition is not ful- 
filled as they are rather simple and give some interesting features of the theory. 

At first we shall determine the energy states of the equation of the meson 
field. This equation is obtained from the Hamitronran in the usual way.” 
The time independent equation is 


(A+h—yQ)U (a, y,2) =0 
8 
where k = | (ho? —x* and FH is the energy of the state. We assume that 
AG, 
the source function is spherically symmetrical and that, if we consider two 
sources, the distance between their centres is zero. Then we can introduce 
spherical polar coordinates 7, 0, p, taking the centre of the source as the 
origin of the coordinate system. Separating the coordinates we have the well- 
known solution 


1 
U =— U; (r) Nihes (0, (). 


5 
U;(r) is determined by 


( a i i y) y 20) Uj (r) = 0. (2) 


dr? y 


1 WENTZEL assumes 4 H = yf U* Od vf U QdV and Pavutt and Hu assume expressions 


corresponding to this one. A ae 

2 Generally the field equation ought to be deduced from the Drrac equation and no 
from the Hamironran in Chapter I, as this HamMIronrANn is usually not correct except to 
the second approximation (cf. § 2 of this chapter). In the particular case here, the Hamiu- 


TONIAN is quite correct. 
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Corresponding to a given value of / we have 2/ + 1 values of m, i.e. the | 
corresponding energy state shows a (2/ + 1): fold degeneracy. Now we imagine 
the whole system included in a sphere with the same centre as the source, | 
requiring U to be finite, continuous and to have a continuous derivative inside | 
the sphere and to be zero on the surface of the sphere. Thus U must be:='08)| 
when r—0. The radius R of the sphere may be considered to be large. | 
Summing the moduli of the energies, we let the radius of the sphere tend to 
infinity. Although the total field energy is a function of this radius, the inter- 
action energy is independent of it. | 

If 2 =0 the solution which is = 0, when r =O, is 


Uy (BR) = e Vig Jizy (ko?) 


where c is a normalization constant and J7+3 (k 97) is an ordinary BrssEL func- 
tion. When r=—k 


er 1 a; 
Us ( hy = (") sin (ky R — 412) +0 (;,) == ((r (3) 
Neglecting terms of the order of magnitude of * this gives 
kR—tla=n2a (4) 


where n is any integer such that kp is positive. We have 
Eox = ai h eV + x2. 


As we are only interested in the modulus we may drop the minus sign. If 
2+0 and if it tends to zero swiftly enough when 7 tends to infinity we 
have approximately BrssEL’s equation for large r, i. e. 

| 


U (R) ~ sin (bR—41x + 5,1, y) + 0(3) 6) 


where the phase shift 6(k,/,y) is determined by the particular choise of Q. 
Neglecting as before, U(R) =0 gives 


kKR—tla+ 6(k, ly) =nx=khR—tiaz (6) 


or 


Dropping the minus sign as above, we get 


Ey =heVR + Pano|/ x L972 2 ky Ho bD) 


e ‘ il 
where we have neglected terms of the order of magnitude of R under the 


radical sign. 
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iL Ko | 
RV a toty) +0( zs)! (8) 


The total self-energy of the meson field is 


Be = hel Vis vi 


foe) 


For= >) (21+ 1) >) Er 


i=0 ko 


where the sum has to be taken over every value of ky determined in (4). 
The total energy of the whole system is 


=) Pei ess (21 + 1) 2 Ei 
i=0 k 1=0 

where & is defined as a function of ky in (7). If we keep J constant we can 
put Ak: 
and the field 

A Bp = Br— Bon =" > (21+ 1) >) (Ex — Ex) Aky = 


1=0 ko 


=1 from (4). Thus the interaction energy between the source 


hex Ko 
serena (21 a= ll 2 Vier. ab (ko, J, y) Ako 
and in the limit of an infinite radius of the sphere 
AE=lim AE es (21 +1) i ho = 5 (ho, 1, y) dk (9) 
2s ie ae - 
R=0 Ea =O VE + x7 ae 4 
We shall discuss three cases: 2 = _ (A), Q=5 (B), and 2 = constant, 


r 
f— a given radius a and Q=0; r>a (C). In ae two first cases, however, 
[aQav is not finite. The solution of (2) when r= R is in (A) not exactly of 
the type (5) and we cannot determine the self-energy in this case. 

We shall determine the self-energy of the source 4 #(y) and the interaction 
energy between two identical, coincident sources. We get this interaction 
energy if we subtract the self-energy (2 4 H(y)) of two sources, both with the 
interaction constant y, considering their ne to be infinite, from the self- 
energy AE (2+) of a source with the interaction constant 2 y — which source 
we get if the two sources coincide. The linear terms in y disappear and there- 
fore we are particularly interested in the higher terms. 


(A) If we put Q= : we get 
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The solution which is = 0 at the origin is given by the confluent hyper- 
geometric function’ M ; 4 (2kir) and for r=R we get 


BE 
: y Ie 
U (R) ~ sin kR— 57 log 2hR—1, + 4(b,1,y) (5 A) 
where 
5(k, 1, y) = arg r(i+1+i%). 
We have 


eee ee? log 2hypR — d(Ko, I, ) | 
° 2ko R R 


In the expressions for the self-energy of the source the term containing log R 
tends to infinity with R. This term, however, is only linear in y and does not. 
give any contribution to the interaction energy. Omitting this term we have (9) 


ig ike es re ee Fe) 
< AE= Ser y | gto we rie 1+ igh diy. 
We can expand? 
~<a ye 
as) aN [ee ity plo |e 
og (+14 is )= > ie hen can tg To Dhatar wee 


Omitting some terms linear in y, by an elementary integration we get 


Bt ie. betsy V ( y ) 1 Sye areal 
A 21+1 pee oe eae aod Sy ee 
ch paras 12 n) OTN enh Oe Sala 


For 1>L and n>WN (£ and WN large) we can expand the terms in the series 
in a Taylor series 


i. : § 
= | y ee ey 1 m | 
Alay (i aa x creel arccos = 'e ie) 


2%ltn 


S yf 2 | Vesa" \s oleae i eae 
eo ron (eg png ia (sa) 4 


but 


‘ See e.g. WHITTAKER-Watson, Modern Analysis 1940 Ch. XVI. 
Nietsen, Theorie der Gammafunktion 1906 p. 23. 
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Thus the double series does not converge. In this case we do not only get 
an infinite self-energy (of which we cannot have any explicit expression) but 
an infinite interaction energy too, the divergency of which is linear. 


(B). If we put Q=% we get 


a +7(04+ 1 
(+8 P um =o. (2 B) 


Similarly to (2) its solutions are Brssex functions the order of which are 
ele t pe i SEAS ey 
Cees erp i+ gt 
From (5) we get 
dL = FU —) 


and from (9) 


ch K=0 j_-9 2 0 ks via 
K 
S = AU 
=— li 21+ 1)=-(’—1 ky + 
K=oo 21 ( 9 ( 8 Vig+ x ° 


° { \ IU ko 
=— lim: 20 De (Ee = 1k te 
ee le ( ) 9 ( Es Vie ee) 0 
K K 
TX 0 lo 1 | ko 
= aa df - =dky + 
Malla oY oi Gis ie! Var ey 


Thus a finite K would give an infinite self-energy but a finite interaction 
energy. When K tends to infinity the interaction energy tends linearly to 
infinity. 

C). We put 

“ Q=1 rsa) 


Q=0 r>as 
and get the following differential equation 
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iy, Siar nag | 
ant (HS OM 


2 
ad Be (ee al ee "oe U(r) =0 valid for r>a 
r 


=( valid for r=a 


(2 C) | 


yr 
where 


ky =VR—y. 


These equations too are transformed BrssEt’s' equations and their solutions — 


are Vr times Brsset functions of the order of 1+ 4. As U must be = 0 at 
the origin, the inner solution (where 7 Sa) is 


U; (r) ~ Viyr “Ji43 (ky 7): 


The outer solution (r >a) is 


U.(r) ~ Vkr [Jita (kr) — ¢ Yin (k7)] 
where 


Ving (hr) = (— 1) J ae?) 


is Wrper’s function and c¢ is a constant that we shall determine. The con- 


tinuity of the solution and its first derivative requires: 


Fe yy nae kere / Ue aN 


where 
dUi(r) _ _hy pe eg ae 
ay "Ven Ldi+4 (kyr) + kyr Silky) 
dU. (r k 
= es ve! L(Jins hn —e Yuko) = £7 Jia) — ea Gon 


From these equations we can derive ¢ and from (5), (3) and 


VER Yiey (kR) ~ ai ae (Zc = 13 
we get 
6 (k,l, y) = arctg ec. 
We thus get the energy (9) 


ioe} 


S ” ky dk 
a =e rD | a8, 


VRE + x? 


oe 3 (Ky a) S143 (Koy ‘ 
koa Wipes + (ko a) S144 (ko 0,4 


d= koran ae 
a) — ko, @ St + (Koy 


* Watson, Theory of BrssrL Functions, l:st ed., 1922. 


a) Ji+y (ko @) 
a) Visa (Ko a) 


> arctg 


(9 C) 
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with 
ko, = Vik — y. 
In order to investigate the convergence of this expression we can divide it 
into four kinds of mtegrals 
K lore) 


=> | => f => i 


0 


and (0) = y if 


(a) is convergent as it has finite limits and its integrand is continuous (at 
least if y is not too large). 
(6) is convergent. In fact? 


J (1 (tgh a—a) (1 o(*)) 
Miko OS eee | 
: V22y tgh a 2 


e—? (tgh a—a) 1 
Ye (kya) ~ a (1 aE o(*)) 
Vinyteha y 


where a = arccosh i and similarly for 
0 
dh (Koy a). 


These asymptotic expansions are valid for large » and aky = v(1 — «) where 
é is any number > 0. 
Introducing these expansions in (9 C) the integrand will have as factor 


e—2l(e—tgh ) 
The terms in the infinite series of (6) thus tend exponentially towards 0 and 
we get a convergent series. We may even let K tend to infinity. 


If ky and ko, are large we can use LoMMEL’s expansions of the BrssEr 
functions? 


(4a 
(a koy)~” Jy (ako) = ga G ae (a ko)~ POI ates (a ko). 


Expanding too the arctgc and Vk3 + #2 we get for (y) and (6) the integrand 
oa AE 


2 
= (1 — Le + | ye (Ji-4 Gh) Jy, 3 (@ ho) — Ji43 (a ky”) + 


she (y a)? le ae (Ji 4 (a ko) Jy 4.5 (a ho) — Tiss (akg) S,, 3a ko) oad 
=) (Ji 3 (@ he) Ty, (@) — Faz y @ ho?) (Vis (hy) Jy, 3 (a ho) — 
Sy Soman ky) | +- |. 


1 Warson, loc. cit. § 8: 4. 2Mlocaicite $95: 22. 
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Assuming in (y) that ky > 1 we can introduce here the expressions of Ji+3 | 


and Y +34 in a finite series and stop at the second term. 


ee eee I 
Ji+s (ako) -/ sin (1% 4 | cos (at —H2) + ve 


a ko 2a ko 


oe La\ To PIs : =) 
Yas (ake) =—|/ = cos (ic =) Vee sin (1% 9 I: 


and we get the integrand = 


_ af, 2 : =|) — Ae yop (s- 1 aie 
~ 9 \r@ A 1 ro(p Cee ea ps 4 — 2 cos” (ako ayy 


We have thus seen that the terms of AZ linear in y have a logarithmic diver- 


gence, while the other terms converge. Thus from (y) we get an infinite self- 


energy but the contribution to the interaction energy is finite. This is true 


even if LZ tends to infinity. 

In (6) we cannot stop the expansion of the BrssreL functions in the same 
way as in (y). It is, however, possible to integrate one of the y?-terms of AH 
We have 

K' 
lim [Fuss ee ae 
K'=00, ‘ k 
IK 


K 
ak td, 43 (Gk) Sota (ak) — Fura (Gk) SF, 2 (ak)} 
ep (2) (ee) 


K 


| 


7 sin (vy — p) Ze 
es ALAA pont 2 
( 


(eee eel rm (u — v)(u +» +1) 
| aK {J .3 (a K) J,+4 (a K) — Furs (aK) I, 3 (a K)} PCHGEE 
(E> Via vk) fe pat j 


where we have introduced the expansion of Ji;; (ak) given above. Thus 


< ‘dk 
2D OEEA [ Fg U8 Ako Fy 5 ao) — Tay ake) Ty, 9 (ao) = 
K 


20 ‘ie, 1 
ee <a Oo) | ? * a l 
a 2 Blan ig 


-+- terms containing Bressex functions 


\ 
J 


* Watson, loc. cit. p. 135. 
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As we saw in (8) the terms containing Besse functions tend exponentially 
towards 0. We get the series 


e? ———— which is not convergent. 


We see that in (5) at least one of the y?-terms diverges. Thus even the 
interaction energy has at least a linear divergency, arising from the region 
where ak and / are roughly of the same order of magnitude and both tend 
to infinity. 

We see that the assumption of a source of finite extension in ordinary 
space does not give either a finite self-energy or a finite interaction energy 
— at least not if the sources coincide. Thus we have to make a more radical 
assumption than this one to get rid of the divergency difficulties in the pair 
theory even if we are not interested in the relativistic invariance of the as- 
sumption or in the convergence of the self-energy problem. We must assume 
that for some reason ak and / are not allowed to tend to infinity, at least 
not at the same time and in the same way, i.e. we must introduce some 
kind of “cutting-off’’. 


§ 2. A criterion for the validity of the pertubation theory 


As above we shall assume fixed, extended sources and in addition some 
kind of “cutting-off’’ rule that makes the total interaction energy converge. 
We shall, however, see that if we restrict ourselves to a perturbation of the 
second order, this rule need not be stated very definitely. We can get a 
criterion for the validity of the perturbation theory if, for a fixed value of 
the unperturbed energy, we calculate the corresponding value of the interaction 
energy using the perturbation of the second order, and compare this value 
with that obtained by including the perturbation of the third order. A cri- 
terion for scalar mesons and scalar interaction has been derived by WENTZEL [42] 
and the corresponding criterion assuming another spin-independent interaction 
has been derived by Pauiti and Hu [45]. The latter authors have derived a 
criterion for vector mesons too, assuming a spin-independent interaction. The 
kind of extended source they assume is, however, not the same as ours (see 
note 1 p. 147). Kuiern? has derived the same criterion as WENTZEL using the 
same kind of extended source as we. In this paper we shall use a modification 
of his method, limiting ourselves to one spin-dependent interaction in the 
case of scalar mesons and one in the case of vector mesons. The criterions 
for the other interactions can be derived in the same way. Starting from the 
generalized Drrac equations for the meson and eliminating the components of 
the field function, with the exception of a scalar in the scalar case and a vector 
in the vector case, we can get the KLEIN-GoRDON equation and the Proca equation 
(with three components only) respectively. The right members of the equations 
are of course not zero any longer, but some function of space, homogeneous in 
the field function. Starting from a set of solutions of the unperturbed KiEry- 


1 Not published. I want to thank Professor Kirin for his kindness of giving me the 
opportunity of seeing this paper. 
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Gorpon and Proca equations, respectively, we can, with the aid of a GREEN 
function, write the solution of the differential equation as a power series of 
the interaction constant. By this method we can obtain the values of the 
interaction energy as a function of the unperturbed energy in any approxima- 
tion. Half the sum of the moduli of these energy values give the total inter- 
action energy (cf. Ch. II § 1 and Hsatmars [49]). As the unperturbed energy 
is degenerated we ought to have used a set of solutions, particularly chosen | 
for the perturbation terms. We shall see, however, that this sum is independent 
of our particular choice of solutions and we have, therefore, a complete freedom 
of choosing any set of unperturbed solutions. 

We put the wave function of the meson field = U (r,t), being a scalar, if 
we have scalar mesons, and a vector, if we have vector mesons. The source 
function is described by y 2, where y is the interaction constant. Q is of 
course a function of space but it may be some kind of operator too, operating 


1.0 aa : 
upon U, as ‘e.g. grad (or 5 Og Of mixing the different components of U. It 


may even form a series in ascending powers of y. 
The field equation is 


—2@U=yQU (12) 


which is a scalar equation (one equation) in the case of scalar mesons and a 
vector equation (three simultaneous equations corresponding to the three dif- 
ferrent components of U) in the case of vector mesons. In this case, usually 
the three components are mixed in the right members of the equation. 

We imagine the system included in a cube whose side is very large = A 
and whose volume A*=V. Then we get the values of the energy states of 
the system (#) if we substitute 


Til pitta (13) 
n?(12): 
Thus the field equation is 
E 2 
(4 + (5) — ) U (r) = y Q(r) U(r) (14) 


where £22 may also be a function of #. 
A solution of the unperturbed field equation is 


i Ep 


? ep — Eo 
Vor» ir t) => Vor» (r) é h => Dorr 2 (1 h ) (15) 


where y= 1, 2,3 corresponds to the three different fundamental solutions, 
corresponding to a single value of k in the vector case. If 


[ Odin (r) Vor» (r) d V can Ow» (15 A) 
Vi 


we have an orthogonal and normalized fundamental system of solutions. 
Putting (15) in (12) with y=0 we get 
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Ey =chV|k[? + 22 (16) 


and thus #) is a function of |k| only and it is degenerated with respect to 
the direction of k and in the vector case with respect to » too. For every 
value of |k| there are two symmetrical values of Ey. As above, we have 
dropped the negative value of Ey. Uox» is a constant scalar or vector which 
we assume to be suitably chosen as regards the perturbation. k is determined 
by the fact that Uox, shall be a periodic function of the side of the cube 
(this is an unessential modification of our previous condition, that the wave 
function must be zero on the surface of the space, permitting us to use 
exponential functions instead of trigonometrical ones). The real solution U;,, (r) 
must of course also be a periodic function of the cube. As we assume that 
the sources are concentrated in a relatively small region around the centre of 
the cube, the difference between U,,(r) and Uox,(r) must be small, if we are 
far away from the centre. 

We multiply (14) from the left by Uox,(r) and integrate over the volume 
of the cube. By partial integrations, bearing in mind that the surface integrals 
are —0O on account of the periodicity of U;,.(r) and Uox,(r), we get, if we 
omit for the moment the indices k and », 


[ Us (r)(4U@) dV = [ (U8) A) U(r) a. (17) 
V V 


At the same time we shall expand HL, U(r) in ascending powers of y and write 
down the possible expansion of {(r) i.e. 


E=£y,+ yh, + YB, + yb, +--: (18): 
U (r) = Uy (r) + y U; (r) +97 U2 (r) ++ (19) 
and 
Or) = Qy(r) + 72, (r) + y? Qs (r) foe, (20), 
Thus we get 


(ch)? 
E, E 
+ y? US (r) lee (2 Ey EF, + Ei) Ug(r) + 2 The? U,(r) — 22 (r) Uo (r) — 
= 24(r) U(r) | + 
1 : 
+ y° Uo (r) [cece Bo Bs + By By) Uo (r) + (ioe 2 Ho Be + Ei) Uy (r) + 


“eran U(r) — 24 (r) Tor) — 2a (r) Ua lr) 24 (7) Wer) | +-\av=0. (21) 


(21) shall be valid for every value of y i.e. the integral of every bracket in 
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it must be =0. As Ué(r) is a solution of the unperturbed equation the first 
bracket is identically —0. Remembering (15 A) and (16) we get the energy 


h 1 . iat ares 
zs via al U6 (r) (Q,(r) Up) d V 
V 


Be 


E\ r 
B,=— 9-H | 0800) U,(r)dV ai 
V 


= —_ | | U6 (r) (2, (7) Ur) dV + [ Us (7) (Qo(r) Uo (r)) d v| (22) 
v Ss Vv 


By =— PBs 2 +2 f g(r U(r —By | U8 (r) U(r) + 
Eo 2 Ey . 2 
V 


V 


2 Vk2 + 32 


} | Us (r) (Q,(r) Ua) dV + [ Us (r) (Qy(r) Uy(r) dV + 
x a3 
a | Uo (r) (Q3(r) Ug (r)) d r| 


4 


As we are only interested in the interaction energy of the sources we may 
neglect the self-energy terms, 1. e. all terms which do not contain the distances 


between the sources such as il Uo (r) (Qy(r) Up) dV. As Uy (r), Uz (r) etc. are 


} 
mainly of the same order of magnitude as U,(r) inside the region of the 
sources and small outside it, some of the integrals containing these functions 
may easily be seen to be negligibly small. 


Neglecting in the energy all terms smaller than we get from (22) the 


i 
We 
interaction energy of the first, second and third order of perturbation 


ch i! a 
BE 5 Vig i Uo (r) (Q1(r) U(r) dV eo 
V 
ch 1 ny " . 
AN pe 5 VB | | vam (2, (r) Ug(r))dV + [ove (2, (r) U(r) iv | (24) 
V 4 
with 


AE=yAEK,+YAE,+yAE+-::. 


The AE:s have, of course, generally the indices k, y too. We obtain the total 
interaction energy in the approximation of a perturbation of the second and 
third order, if we take the sums of A Fox, and A Espy, summing for every 
value of k and », remembering our cutting-off rule, and letting the side of 


the cube tend to infinity. In (23) and (24) some self-energies are still left, 
which we shall neglect. 
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In order to determine U (r) we put the expansions (18), (19), and (20) into (14) 
AG Vat 


‘(Uolr) + yUy or) + YUL) +-:-) =0 (25) 


(4 + (2) 3) oro 


which is satisfied by (15) and (16) and 
(eee (ae AB 
(2 (8) —) 0s = (a9 2B) 
Hy + 2H, Ey 


fe (2, (r) (ch)? U(r). 


Neglecting in the right members the terms that disappear when V tends to 
infinity we get 


or 


(A + k*) Uy (r) = Q (r) U(r) 


(A + BY) Uy (r) = Q,(r) Uy (r) + Qp(r) Uo (r). eo) 


A solution of the equation 
(4 + h*) U(r) =f (r) 
where the function /(r) is known, is 
2 [ Gr (r,r’) f(r’) dV". 
If we put the Green function 


ein" (r—r') 
Gy (r,r’) = a Me (27) 


the condition that U(r) shall be penodien in the cube is fulfilled. Remembering 
that 


a Mea g(r) dV! = f(r) 


we can easily see that we have got a solution of the differential equation. 
Letting the cube tend to infinity, we have 


co +1 2 


; il ce fe ” ik’ Jr—r' |w ; 
G (rr) =lim Grier) = 9 a | | | qa par Wd dudg 
oO” —1'% 0 


where r—r’ is the polar axis and u = cos #. 
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dk’. 


Grr) = hk? — k’? 


ovo} 
1 1 (ieee 
2 x | r—r | é 
0 
Taking the principal value of this integral, we get 


| 


(28) 


1 cosk |r — ra 


4n |r—r'| 


We see that if f(r) +0 only within a limited region lim U (r) = 0 as we have 


required above. Thus we get from (26) 


U,(r) = { Gy (r,r’) (Q) 7) Ug’) dV’ (29) 


and 


Us (r) = f Gy (r,r')(Qy(r) U(r) dV’ + 
V 


+ [Gr (rr) (Qq0r') U(r) dV’ = (30) 


V 


=f [ Grr) [Q(r') Gr (rr) (Qi) Up r)] dV dV" + 
VV 

te ip Gri(r,r)(2Q,7) 0507) dV; 
V 


The criterion mentioned above is obtained from the requirement that 
A ber, > y A Esny (31) 


where A fox, is given by (23, 29) and A H#3,, by (24, 30). As the second 
term of (24) is analogous to the part of the first term of (24) arising from the 
second term of U,(r) in (30), we may omit this term when we determine the 
criterion. 

We shall now show that it is permissible to choose any fundamental set of 
solutions of the unperturbed equations if we want to calculate the total inter- 
action energy. We shall here limit ourselves to the perturbation of the se- 
cond order, but the proof can easily be generalized to any perturbation of 
higher order. 

Uoi(r) is an arbitrary ortho-normalized fundamental set of solutions of the 
unperturbed Kiern-GorpoN or Proca equation and Uoi(r) is another ortho- 
normalized fundamental set of solutions which is particularly chosen for the 
perturbation calculation. Both sets correspond to the same energy Eo x |: 

Thus we can put 


Ualr) = > a2 Voi (r) (32) 


where cj; are constants. 
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The ortho-normalization condition gives 
J Von (r )Ua(r)dV = > cme cu | Ub r )Uoi(r)dV = 


ua 
— = Cru C1, 6 ti = cae Craik Cha = Omi (33) 


wh 
and correspondingly 


De Cae Cj, — Ont (34) 
m 


From (23) and (29) we get for the total interaction energy 


he 1 
Ato — >, A Ea, = > , 
ik mu 2 VE + 2 


ff War (r’) [Qs (r') @ (r', r”) (Qi) Vale’) dV dV" (35) 


and from (32) 


: he I! 
A Eo; = Cin Ol aS ae 
ih 2a 2 VR + 92 


‘ff Usi(r ‘) (2Q.(r G(r % r’’) (Q, (r’’) Uou(r’’) | dV’ d Wie 
From (34) we get 


he i 
ery 
ras 2 Ve + 2 


ff Ub.(r’)[Qy(r VE, r”) (Qe) Toile’ dV' dV". (36) 


Thus we are allowed to sum, assuming an arbitrary fundamental set of solu- 
tions of the unperturbed equation, which of course can always be made ortho- 
gonal and normalized. 

With this method we shall now deduce a criterion for the validity of the 
perturbation theory in some particular cases. 


A. Scalar Mesons with spin-dependent interaction 


If we assume the interaction term to be built up of antisymmetric tensors, 
the generalized Dirac equation of scalar mesons is (cf. Hsaumars [49] 110, 
192 and 331) 


(> Dr Pe + vpe— inne Sonny hes Tle 0. (37) 


k=1 
k+l 


1, 2,3, 4,5 


We shall assume the source to be at rest and thus neglect the smaller terms 
of the source function, i.e. the sums in the interaction terms may be taken 
between 1 and 3 instead of between 1 and 4. 


161 


OLOF BRULIN, Meson pair theory | 


We have further 
VYVEVY = VY vavevip =—t8s (k= 1=1, 2, 3) (ef. Ch. I (3)). (38) 


The element of the matrix f; is (cf. Hsatmars [49] equ. 80) 


(m | Bu| 2) = 4 (Om On — Sms Orn). (39) | 
Introducing this in (37) we get 
W112 Obes Ue Sean r1 = 0 
or ee 

pU + weU,= type (Ss X U,) (40) 


where Us is put =U and the vector U, has the components U; (1 = 1, 2 
and 3). 
Further we get from (37) 


p,U + ucU,=0 (41) 
and 


p: Ga 7,0,— pel =], (42) 


From (40) we can express U, in U. As y is assumed to be small we shall 
do this as a power series in y. By iteration we get 


Up U NSU 
1 7 
Tae pe LPS AC SBE) eae ae (43) 


Introducing the series (43) and the expression for U, obtained from (41) into 
(42) we get 
PlpU —ipcy(S3XU,)] + piU + wWekU =0 


or with p=— ch grad, m=—* = 


(4 U-s sa —*) ig 
=—xy div (Sx U,) =1y div (S3x grad U) +? div (S;(S3X grad U)). (44) | 
We have thus written the field equation in the form of equation (12) with 
Q, = 1 div (S3 X grad ) 
Q2 = div (S3 X (S3 X grad )). (45) 


We shall first give the value of the interaction-energy in the approximation 
of the second order. We have already calculated this value in Chapter I and | 
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_ we shall show how to obtain the expression given there. Our assumptions 


regarding the source function are the same as in Chapter I. Thus we get ((15), 
(20), (23), (27), (29), (45)) 


ch 1 ald ioten 


ihe!’ (r'—9!’) 
aout e ld. 0 div’ (s, (r Rex grad’ ae seal 


+4 div” (Ss ér’’) x grad” Wer ern) dvV'd Vv" 


where grad’ etc. means derivation with respect to a’ etc. 


By partial integration we get 


PEST ae | 

; A bij Sl SS v= —— 5 Vow’ Vox pa p98 p2 
i gradee3 ©?” yp (r') grade” | a VY" 
V 


J grad ere (Sa(r") orad &® | dV". (46) 


If Uox(r) is ortho-normalized, we may — as we have seen — chose Uox quite 


arbitrarily. We shall put 


Uon = (47) 


1 
We see that, except for the factor —————., 
Vik? + 2 


sign if we exchange k’ and k’” and at the same time r’ and r”’. Thus we 


every term in (46) changes its 


' can write every term in (46) as a difference between two terms both half as 


large as the original one, having performed this exchange in the last term. If 


we extend the first term by Vk’? + +2 and the last term by Vk? + x7, we 
can perform the subtraction and abbreviate with Vk’? + #2 —Vk? + x. 

Carrying out the differentiations in (46) and rearranging some of the terms 
in the scalar triple products we get 


ABs = 4S V' dk' dk’ - 
ei (k'—k"’) Qr''—r’) SS ee = i! ee ee, 
Vk + ®VE? + 9? (Vie? + x? + VI? + 2?) 
Bt S3(r’)-(h’ X k’)} (S3(r") (kX k’+)}. (48) 
_If we put 
=—8 G. fhe. (49) 
aoe ee: 


vy? AHz; is easily seen to be identical with S33’, in Chapter I. We can for 


: lo carry out the integrations with respect to y’ and gy”. Then we get 
(68 S33") of Chapter I. 
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We saw in Chapter I, which cutting-off rule we had to introduce in order 
to get convergent integrals in (48). If, however, the source functions defined — 
later on in this chapter cut each others, this rule is not sufficient any longer as_ 
the integrals obtained do not converge in this case. In any case, if the dis- 
tance between the sources is very small, the whole theory can not be assumed _ 
to be very satisfactory. Thus in our numerical applications we shall have to 
make an a posteriori assumption concerning the interaction energy inside | 
some small distance, the ‘‘cutting-off distance’, which is of course of the same 
order of magnitude as the “radius’’ of the source defined later on in this chapter. 

We shall now investigate the criterion for the validity of the perturbation 
theory. Owing to the cutting-off for large k, we only need to determine the 
energy for a fixed value of k which we may assume to be rather small. At 
present, however, we have to make some explicit assumptions regarding the 
source. We shall assume that the spin of the source is constant and that we 
have a spherically symmetrical source function (“‘radius’” =a). Thus we put 


| 
| 
| 
| 
| 


S3(r)= out (lr—rul), ff dr—ral)dV=1 (50) 
V 


-& 


where the constant vector r, determines the position of the w:th source anc 
ia (lr—r,|)=0 when |r—r,| >a. 

For the sake of simplicity we shall put f.(/r—ry|) = 6(|r—ryl) so far as 
we get a convergent result in that way. In (23) and (24) we may thus by (45) 
introduce a 6-function and evaluate the integrals. From (31) it follows that 
the criterion is given by 


Vir (rx) > Y Vox (r,.) (51) 


where we may particularly put r,—r,. In our calculations we shall assume 
the distance between the sources to be >a. We shall only consider the case 
of two sources S, and Sz. 

From (28), (29), (30) and (45), we get (omitting the index k in Ux(r)) 


U; (r,) == <_s div’ (Ss (r’) x grad’ Us (r’)) d Kite 


[acme tee 


1 | COS RIT 1 
ree 


If we omit the terms that do not contain |r, —r,| we get with the aid of | 
(15) and (50) 


1 cos k|r,—r| 
A U = = 1 Ler 
1 (71) re ee bryce Pee x ke'*"2 Vor. (52) | 
As k is small we get the dominating term proportional to rae ral Further: 
| ics ae a 
we get ‘ ‘ 
: l ale cos k|r, —r’| 
U2 r,) =— — + | | if yw Lge sant ’ 
(r1) (i eet adV'dy eee div’ | S,(r’) x 


cosh | rir | 


Xx grad’ foeerin div’ (S3(r’’) X grad o'r | (53) ) 


where U3(r,) is the first term in (30). 
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Terms where r’ and r’ are both in S,, do not contain ry i.e. they may be 
omitted, as they do not contain the distance between the sources. Terms where 


t 


| r’ is in S; and r’” in S, or vice versa, may be neglected as they contain higher 


i. 1 : : : 

| powers than the second of ered for if we insert in these terms the values 
the he 

of S;(r’) and S3(r’’) respectively, being essentially 5-functions, we get through 


=e : 1 Si d. 
' partial integration factors as grad --—~——,, multiplied by grad ———— where 
=, | | he r,| =P | Sup | 
a ih aie a 4 
Fr is in S,, or grad’’- grad” reer where r’ is in S;. Thus we may limit 
Mah = 5 
ourselves to the terms where both r’ and r’ are in Sj. Putting S3(r’) as a 
6-function and carrying out the integration with regard to r’ we have 


re iow 


AU) =a | dv” grad cos k|r,;—r fo, s 
iS 


et I 


r=Pr, |r—r’| i 


From (50) we put S3(r”’) =f (r) 62, r=|r” — rg. 
Further we put 


div (S3(r”) X ke” Uox). (54) 


a k|ry—r'| 
(42)? van [ti — 7" | 


Oo xk etkre Vox eee 4 


X 02 =, 


Neglecting the variation of the factors containing k inside the source we get 


AUe(r;) = i| dV (c - grad ~) div (f(r) w) = 


ifay (v- rad 5) (w - grad f(r) = (55) 
=i f aves w-Ff (7) 


Putting for example the polar axis parallel to v we get 
, 47 fe : An 
A Ua iete tee i (r)dr=1 3 wv (0). 
0 


If we assume that the source has constant density and the radius a i.e. if 


we put 
f(r) =4O=—3 (56) 


47a 
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we get 


A Us (r;) eat SO (0. x grad cos bl rites “l) (6, x k e'#2), (57) 


(4 mt)" a r=1 | Lisi om. r| 


As o is of the order of magnitude of 1 the criterion for the validity of the | 


perturbation calculation runs 


pee 
40° 
or from (49) Sane 
4 83 é€ ~ 
LSE” ae 9 
x a® hie : (59) 


This is the same criterion as Pautt and Hv obtained. 
Inserting (28) and (45) in the second term of (30) we get 


1 Oe | , Lvs a { So (r’)X (Sg (r’) X grad Uo(r’))} dV’. (60) 
4x lr; — r'| 


As above the integration is carried out only over Sy. One of the factors S3 (r’) 
can be put equal to 626(r’—rg) and the other = 63/(|r’ — rg|). Carrying 
out the integration of (60) we get 


A U2 (r;) — ‘ grad ear Nias ae Sra -£ (0) (Gg x ‘09 x k) Uone*”2. (61) 


re a ea PRE 


3 
Putting / (0) = ~—- and comparing (61) and (52) we get 


4a 
3 
ek 
or from (51) 
Os3 é 
| a? he “ (62) 


which is about the same criterion as (59). 


B. Vector mesons with spin-dependent interaction 


We assume also now that the interaction term is built up of antisymmetric 
tensors. The calculation is mainly the same as in the scalar case. The Drrac 
equation runs (HsaLmars [49] 167, 192, 363) 


4 4 
{> re Br + iwe—iype > PpenyBr Bros! v| aw (63) 
k=1 ues J mn 


where (Hsatmars [49] 161) 
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| 
<1 (58) | 
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and (Hsatmars [49] 139. The index combinations 135 shall be used.) 
(mn| Bul m' n’) = (64) 

=4 (Omk Om's Onn! ae Ons Om'k Onn' ae Onk Om's Omn' ar Ons Om'k Oman’ a 


re Omk On’5 Onm' =e Oms On'k Onm’ =i Onk On's Omm! a= Ons On'k: Omm’). 


: We make the same assumptions as regards the source as in (38). Introducing 
| (64) in (63) we get the following equations 


= 7d Ts » Pm Usn — Pn Usm + UC Umn = 0 (65) 
3 
>) pe Via — we Uss = 0 (66) 
k=1 

and 

4 

1, 2,3 >) pe Uin — w ¢Usn + iy we(S3 X U)n = 0 (67) 
| k=1 
| where 
(UC), = Us5n. 


We multiply (66) and (67) by py and p, respectively and sum, remembering 
4 

| that ps Dn Pr Urn = 0 as Uzn is an antisymmetric tensor. We get then 

nm, k=1 


4 3 


>) Pe Use — ee pr(S3 X Ue =0. (68) 


k=1 


Inserting Umn from (65) into (67) we get 
4 4 
i — i 7p 3 ~ Dk Usm = Dn ey Pk Usk se we Ce Usn a ty wc (S3 x U)n = 0. (69) 
k=1 k=1 


With the aid of (68) and introducing the definitions of p; in (44) we get 


2 
(4 — 2s ga — #4) U = ty (ead div — #4) (5, x U) (70) 
Cc 


- Thus we have got the field equation in the form of (12) with 
Q, =1(grad div — x?) (S3 X ) (71) 
2, = 0. 

As above we get the interaction energy (from 23, 27, 29, 71) 
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: ‘ Det sey ee Mare 
| / Ucipen "se | (eraa div’ — x?) (s, (r)x% LESS oa) 


-(grad”” div” — x”) (S3(r") X Uor’s oun | dVdV (72 | 


In order to evaluate these integrals we shall use the following vector relations. | 
By partial integration and assuming that all surface integrals vanish we have! | 


fa-f( r) grad div F(r)dV = 
= — fa: grad f(r) div F(r)dV = 
= | (a grad - grad f(r) - F(r)dV 


and 


[ f(r) grad div F(r)dV = 
= — grad f(r) div F(r)dV = 
= [Fo grad - grad f(r)dV. 


As S;(r) = 0 outside S, and S, these formulas may be used. We may put (cf. 36) 


3 
0k3 = ae 


Uor1 = alee Uox2 = Va VV 


Yh 
VE Vy 


assuming the same coordinate system as in Chapter I. 
Using first (73) and then (74) we get 


he 1 1 
A Hs; = lim ye ms > 9 Vie ae 


Keka 


: sea aV'dV" |x, (grad grad e-"" = x7 e-1©* )]- [Sadr ix 
vv 


x ef" {85 (r'"") X a, er) (grad - grad e780" — 42 e-ik' 1) 


Analogous to (48) we get 


ead i ; ai it 
A Eo; = uci u iL ae Bao eee 
Ex=—- sar | oF ak dV dV" a 


{S83 (r’) X [(a-k) bk’ + eat 


-{[(S3(r") X a) KU] + 22 (Sg (nr) X &)}. 


OF OF OF 


"0am 75 a “2 On" 


1 


a grad: Fir)=a 
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(74) 


(73) 


? 


(75) 


(76) 


(77) 
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If we put 
ee all 


y= baz 3 Cw, (78) 


we get the same value for the interaction energy as in Chapter I V3"3". This 

can easily be seen. We may for example, having introduced the same polar 

coordinates as in Chapter I and exchanged the coordinates as in (48), carry 

out the integrations with respect to y’ and wy’. We get then the same expression 

as in Chapter I (51 V3-3), where we have neglected the terms containing Sy (r). 
We shall now investigate the criterion in the case of vector mesons. We shall 

make the same assumptions as in the scalar case but we may now put k = 0. 
From (15), (28), (29), (30), and (71) we get 


U,(r,;) = (3) firs | (grad’ div’ — x”) (S3(r’) X Uy) dV’ (79) 


and 


SSS (2) i ema {aad de) Is, (eX presto 


-(grad” div” — x2) (Sir) x vu»)! dv'dV" (80) 


As in (52) and using (74) we get 


AU, (ry) = — ra x U5) (saa - grad : he : : (81) 


PosACi Th |r; —r2| 


In U,(r,) we may as above omit terms where r’ and r” are both in S, as 

these terms do not contain the distance between the sources. Further we may 

neglect those terms where r’ is in Sj and r” in S, for if we put S3(r) as a 

6-function and use (74) we see that these terms contain grad - grad ree 

r=P2 i: aon 

b ———— 

? [ri — re | 

Thus we may limit ourselves to discuss the terms where r is in S, or in S, 
and r’ is in Sp. ' 

In the first case we may put S3(r’’)=6,6(r—r.) and S; (r == 644 (2) 


multiplied by itself and thus contain higher powers than the third of 


where r=|r —r,|; in the second case we may put S3(r) = 6,6(r — re) 
and S,(r’) = 62f(r) where r= |r” — rol. 
1 
Neglecting terms of higher power of ie . | than the third or x* times the 
a Te 


first and using (74) we get 


2 
: holt , alee 
AU, (r,;) = — (7) fav , (grad div’ — x?) f(r) - 


eee 
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(A) far’ Mgr bites ee 

——_ — . 1G ~ — —— . 

4a} , a Pe reel |r; —re| 
S 


ie 4 ; {(grad”” div’ — x®) f(r) (6g X Uy) | . 


We have, however, 


grad div (f(r) v) = v grad - grad f(r) 
and 


| ~ grad -grad f(r) dV = 


= [2 saa’ ray + fare eraderd) = 


,, 
| Pere “fi je ie 
E [ 2 Fall r ray Fea he ls 
(if we use polar coordinates with v as the polar axis) 
_ 4x ee i) 2 fre _ 
ae v| = ag fe oe cae LAr bt = 
0 
4 ; r 
=F [tion Ze [im +: anv | f(0) = 
0 0 0 
4a 
Pian te Ol 
3S 
(i we assume /(r) constant = 4 
4a 
Dv 
Pi oat 
and 
[,foav=s2frsma <r 
0 
Thus 


bay 3 20” 
AU, (r,) = (7. ie | al 


il 2 
[x (ioe x Op (ead grad eee eal te 
r=, | Py —_— 179 
+ 40s OOK Oar< U,)} (srad «grad / = eae oi 
Li Ni 1 FTPs 
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Thus we get from (86) and (81) 


aS. SREP Vi 
4na° Paes | oe ine 
or from (78) 
2 Cpe | 28 
x a? (; ) Pa (en) 
and 
Pcs (=) Ss 
| xa \ie} | — 2 = 


If (xa)<V2 (87) is the stronger criterion, which is usually the case. (87) 
corresponds to the criterion in the scalar case with spin-dependent interaction 
and to the criterion found by Pavitt und Hu [45] with a spin-independent inter- 
action in the vector case and (88) corresponds to the criterion that Wentzet [42] 
and Kiri found in the scalar case with scalar interaction. As we saw above 
(87) arises from the terms containing two differentiations with respect to space 
and (88) arises from the terms that contain x” as a factor, i.e. that do not 
contain any differentiations with respect to space. 
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